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Abstract 

A nonrelativistic Bose gas is represented as a grand-canonical ensemble of 
fluctuating closed spacetime strings of arbitrary shape and length. The loops are 
characterized by their string tension and the number of times they wind around the 
imaginary time axis. At the temperature where Bose-Einstein condensation sets 
in, the string tension, being determined by the chemical potential, vanishes, the 
system becomes critical, and the strings proliferate. A comparison with Feynman's 
description in terms of rings of cyclicly permuted bosons shows that the winding 
number of a loop corresponds to the number of particles contained in a ring. 
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The first experimental realization of Bose-Einstein condensation (BEC) in a cloud 
of weakly interacting 87 Rb atoms trapped in a magnetic field by a group in Boulder, 
Colorado and the second definite observation of BEC in a system of Na atoms at 
MIT [^|], turned this field into one of the most active research areas in contemporary 
condensed matter physics. This recent burst of interest motivated us to consider once 
more the fundamental properties of BEC. 

At present there exist at least three different descriptions of this phenomenon. In 
the first, and intuitively probably the simplest one, BEC is understood as the condensa- 
tion of a finite fraction of the total number of particles present in the zero-momentum 
ground state [B]. The second description, due to Bogoliubov [Q], uses the methods of 
quantum field theory [Sh. The formation of a Bose-Einstein condensate is here signaled 
by a nonzero expectation value of the field describing the particles in the system. The 
third description, due to Feynman [Q], is based on path integrals [Q]. BEC is now un- 
derstood as the occurrence of large rings containing many bosons which are cyclicly 
permuted in imaginary time Using a mapping onto a certain problem of clas- 

sical statistical mechanics [|l0] to numerically evaluate the path integral, Ceperley and 
Pollock [ [TTj ] showed that this picture also provides a powerful computational tool with 
which even a strongly interacting system like superfluid 4 He can be accurately de- 
scribed (for a review, see Ref. [|l2|]). The classical system is a chain of beads connected 
to springs. 

In this Letter we discuss an alternative, string picture of BEC. This description is, 
like Feynman's, based on the spacetime approach to quantum mechanics, i.e., on path 
integrals. The specific path-integral representation we arrive at was first put forward 
by Montroll and Ward [|l3|]. Physically, it describes a grand-canonical ensemble of 
fluctuating closed spacetime strings of arbitrary shape and length. We study the critical 
behavior of this loop gas close to the temperature T c where Bose-Einstein condensation 
sets in, and show that the loops proliferate below T c . That is to say, the loop gas 
undergoes a Hagedorn transition at T c . 

The critical behavior we found here is similar to that of vortex loops in superflu- 
ids [ |l4| ] and in type-II superconductors jl^]. (For a detailed account of the so-called 
dual approach to superfluidity and superconductivity in which the vortex loops are the 
elementary excitations, the reader is referred to the textbook by Kleinert [|l6|].) 

Our starting point is the Lagrangian of a nonrelativistic free scalar field theory in 
Euclidean spacetime x = (t, x), with r the imaginary time [|[|: 

£ = q>* Ud T -n- ^V 2 ] <p. (1) 
\ 2m / 

Here, <fi is the field describing the scalar particles, d T — d/dr, [i is the chemical poten- 
tial which accounts for a finite particle number density, and —h 2 V 2 /2m, with m the 
mass of the particles, is the kinetic energy operator. Depending on whether the parti- 
cles are bosons or fermions, the field <f> is an ordinary commuting or an anticommuting 
Grassman field. 

We adopt the imaginary-time approach to thermal field theory, where the time vari- 
able t is taken to be of finite extent < t < fi(3, with (3 the inverse temperature 
j3 = 1/fesT (fee is Bolzmann's constant) [|l7|]. As a result, the energy variable con- 
jugate to t takes on discrete values. Depending on whether <j> describes bosons or 
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fermions, the field is periodic or antiperiodic in r, 4>(hf3, x) = ±0(0, x), and can be 
expanded in a Fourier series as: 

n— — oo v ' 

where k = (w n , k), k ■ x — uo n T — k ■ x, and uj n are the discrete Matsubara frequencies, 

J 2mr/ti(3 for bosons 

(2n+l)7r/n/3 for fermions, ( } 

with n an integer. 

At finite temperature and density, the thermodynamic properties of the system are 
specified by_ the partition function Z. Expressed as a functional integral over the field 
6 it reads 
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up 



D<j)*D<f> exp y-jrj dry d'xCj . (4) 

Since this functional integral is Gaussian, it is easily evaluated to yield 

ln(Z) = In {DctT'' [(-ihuo n - fj, + h 2 k 2 /2m) /%] } 

= -r;Trln [{-ihuj n - fj, + h 2 k 2 /2m) /%] , (5) 

where we used the identity Det(A) = exp[Tr ln(A)]. In Eq. (|j), r\ — ±1 for bosons and 
fermions, respectively, and the trace Tr stands for the integral over Euclidean space- 
time x as well as the summation over Matsubara frequencies uo n and the integral over 
momentum k. 

We wish to derive an alternative representation of the system and express Z not as 
a functional integral over the field (j>, but as a path integral over spacetime loops. In 
the string picture that emerges, the system is characterized by the string tension of the 
loops and the number of times they wind around the imaginary time axis. 

Before considering the partition function, it is expedient to first study the propaga- 
tor G(xi, X2) = G(xi — X2) of the theory, 

G(x) = — T f— - e-*'* (6) 

and write it as a path integral over spacetime trajectories. Employing Schwinger's 
proper-time method [^8[], which is based on Euler's form 

1 1 f°° ds z sa 

-s z e~ sa , (7) 



a z T(z) J 

with r(z) the Gamma function, we write the right-hand side of Eq. (^) as an integral 
over proper time s: 

f°° 1 r ffih 

G{x) = J dse^ h ^J2 e ~^ n{T ^ S) J 7j^3 elkX ^ k2s/2m - ( 8 ) 
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We next use Poisson's summation formula, 



E 



E 5 ( a ~ w ) 



(9) 



to replace the summation over the Matsubara frequencies by another one 

oo 

-i Un (r-,) = g r} w 5[s-(T + w?il3)], 



1 

hp 



oo 

E • 

l — — oo 



(10) 



The physical meaning of the new summation index will become clear shortly. The 
delta function in Eq. ( |To| ) links Schwinger's proper time s to the Euclidean time r. The 
propagator now becomes: 



G(x) 



w=0 



rj exp 



^(r + whp) 
n 



2nh{T + wh/3) 



3/2 



exp 



mx 



2H(t + wHP) 



(11) 

According to the path-integral approach to quantum mechanics [[7]], this can be written 
as a sum over all possible spacetime trajectories running from x(0) = to x(r + 
wh/3) = x: 



G(x) 



w=Q 



with 



G w (x) 



x.(r-\-whf3)—x 



(0)=0 



Dx(r') exp 



dr' L 



(12) 



(13) 



Here, L is the (Euclidean) Lagrangian of a nonrelativistic point particle with mass m 
in a constant potential — [i: 

L = \m± 2 {T)-n, (14) 

where x(r) = <9x/<9t. Remembering the (anti)periodicity in r with period 7i/3, we see 
that the summation index w, which was introduced in Eq. ( 10) to replace the summation 
over the Matsubara frequencies, denotes the number of times a given trajectory winds 
around the (Euclidean) time axis. The sum over the winding number w in Eq. ( p"2| ) then 
implies that the trajectories, which started at x(0) = 0, can wrap arbitrarily many times 
around the time axis before reaching their end point x(t + wh(3) — x. 

We next turn to the partition function. The two main differences with the propagator 
are that, first, instead of a factor 1/A in Eq. (^|) we have a ln(A) in Eq. (Q). And second, 
the expression for the partition function contains, in contrast to that for the propagator, 
no IE-dependence. Taking these differences into account, we obtain 



ln(Z) = r)V 



ds 



d J fc 
(2^ 



iKk. 2 s/2r, 



where we used Euler's form (|7|) in the limit of small z, 

Ha) = - f°— e~ sa 



(15) 



(16) 
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ignoring an irrelevant additive constant. The integration over spacetime produced a 
factor U/3V, with V the volume of the system. Following the steps leading to the 
explicit expression (|ll]) for the propagator, we arrive at the familiar fugacity exp(/3/i) 
series [0] 



w—l 



with 



A = h^/2n(3/m 



(17) 



(18) 



the de Broglie thermal wavelength. In Eq. (|l7|), we omitted the singular zero-tempera- 
ture contribution corresponding to zero winding number w = 0. 

The corresponding path-integral representation of ln(Z) is obtained from Eq. ( fl7| ) 
in a similar way as the path-integral representation ( 12 ) of the propagator was obtained 
from Eq. (JTT[). We find: 



\n(Z) = vY,V w+1 G w (0)/w, 



(19) 



w—l 



where G w was expressed as a path integral in Eq. (O). The extra factor 1/w stems from 
the factor 1/s in the Schwinger representation ( ^q ) of the logarithm, while the zero 
argument in G w (0) reflects the absence of any spacetime dependence in the partition 
function. Explicitly, 



G w (0) 



x(whf3)=0 i ^ 

Dx(r') exp — — 

x(0)=0 V 



wh/3 



dr'L 



(20) 



implying that this path integral is a sum over all possible closed spacetime trajecto- 
ries, or loops, that return to their starting point after a time wh/3. Since one can start 
traversing a closed trajectory anywhere along the loop, the extra factor 1/w (or, more 
precisely, 1 / wfifj) in Eq. ( p^ ) can be intuitively understood as arising to prevent over- 
counting. The arbitrariness of the space coordinate of the starting point yields the 
volume factor V . 

In summary, we expressed the partition function as a path integral over spacetime 
loops. Each loop in the loop gas can be of arbitrary shape and length, and is charac- 
terized by a winding number w, showing how often it wraps around the time axis. The 
sum over w in Eq. ( p^ ) implies that, in principle, this can be arbitrary many times. 

With this insight, we return to the explicit expression (|l7|) of the partition function, 
restricting ourselves to bosons. To be able to carry out numerical studies, we put the 
system in a box with sides of length D, so that V = D 3 . As a result, the integration 
over momentum is replaced by a summation over the discrete variable k„ = (27r/D)n, 
where n denotes the integers (n x , n y ,n z ). This leads in turn to the following change 
in the propagator (111]): 



exp 



2%{t + wh/3) 



E 



exp 



2h(T + whp) 



(x - nD f 



(21) 
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Instead of Eq. (|17|), we now obtain as explicit form for the partition function 
where the prime on the sum will be explained below and 

oo 

w(x)= e ~ xn2 - ( 23 ) 

n— — oo 

The average particle number density n = (N) /V in the box is obtained by differenti- 
ating ln(Z) with respect to the chemical potential: 

In the following we shall assume this density to be fixed. As numerical value we take 
n = 22.22 nm~ 3 , which is the density of liquid 4 He at saturated vapor pressure. When 
varying the size of the system, we always adjust the particle number (N) = N in such a 
way that the density remains constant. For a given volume, Eq. (^) then determines the 
chemical potential as function of N and (3. Figure [l] shows the temperature dependence 
of \x at fixed particle number TV = 10 7 . We observe that it exhibits critical behavior, 
tending to zero as 

M (T)~-(T-T C ) C , (25) 

with the critical exponent c = 1.996. Here, T c is the condensation temperature of a 
free Bose gas, 

2/3 

(26) 



Tr 



2irh 2 



C(3/2) 



where ( is the Riemann zeta-function, with C(3/2) = 2.61238 • • •. Numerically, taking 
the 4 He value also for the mass parameter, m — 6.695 10~ 27 kg, one finds T c = 3.145 
K. 

The exponent c is related to the correlation length exponent v and the dynamic 
exponent z via 

c = 2vz. (27) 

The dynamic exponent appears here because the strings are embedded in spacetime 
and parameterized by Euclidean time. (For strings embedded in space, we would have 
c = 2v instead.) The value c w 2 we found numerically agrees with the expected 
Gaussian exponents v = 1/2 and z = 2. 

Since ln(Z) is written in Eq. (|2|) as a sum over the winding number w, also N can 
be expressed in this way: 

N = Y'(N W ), (28) 



where 
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Figure 1: Fit of the chemical potential. 



denotes the average number of particles contained in loops with winding number w. 
To obtain a loop interpretation of the total particle number N, note that (N w ) is simply 
related to the average number (L w ) of loops with winding number w via 



so that 



(30) 



(31) 



This means that the particle number A*" is given by the total loop "length" measured 
in units of the characteristic "length" scale a — %j3. We used quotation marks here 
because the loops are embedded in spacetime and parameterized by Euclidean time, so 
that a has the dimension of time, not of length. 

Since we keep the number of particles in the box fixed, the summations in Eqs. ( |28| ) 
and d3l| ) have the upper bound N\ there are no particles available to fill longer loops. 
To indicate this we decorated the Sigma's with a prime. 

The partition function written in terms of loops reads 



(32) 



It follows that each loop gives the same contribution to the pressure P = ln(Z)/pV, 
and also to the energy density 



d(f3P) 



d(3 



(33) 



of the system, independent of the loop length or, equivalently, the number of particles 
contained in the loop. 

In the limit of large D, the loop length distribution (BOh with Eq. (E9) reduces to 



m/3\ 3/ V-D x " 
~2n 



-5/2 



-f3al 



(34) 
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Figure 2: Loop length distribution for various temperatures. 



where I = wa is the "length" of the loop measured in units of a — h{3, and a = —fi/h/3 
is the string "tension". Equation U34 ) is the typical form for the loop length distribution 
of Brownian strings fll^]. Figure^ shows this distribution for various temperatures, 
assuming again liquid 4 He values for the mass parameter m — 6.695 10~ 27 kg and 
the (fixed) particle number density n = 22.22 nm -3 . For T > T c , the chemical 
potential is negative and the strings have a positive tension. In this temperature regime, 
we see an exponential decay. When the condensation temperature T c is approached 
from above, the chemical potential and thus the string tension vanishes. The decay 
is seen to change from exponential to algebraic. Because their tension vanishes, the 
strings proliferate. (The temperature where this happens is known in this context as the 
Hagedorn temperature.) BEC, therefore, corresponds to the appearance of long loops, 
wrapping arbitrarily many times around the time axis [^|. 

A similar proliferation of line-like objects we found here, occurs in superfluids and 
type-II superconductors [16|. The line-like objects in these systems are vortices which 
are the elementary excitations of the dual approach. Inparticular, the string tensions 
and the loop length distributions behave similarly [ |l4[ O ]. Only the numerical values 
of the critical exponents differ as a superfluid and a type-II superconductor belong to 
the XY universality class, whereas a free Bose gas has Gaussian exponents. Another 
difference is that a dual transformation interchanges the low- and high-temperature 
phases, so that the vortices proliferate in the high-temperature phase. 

To study numerically the correspondence between BEC and the appearance of long 
or, for an infinite system, infinite loops, we follow Ref. [|l4|| and distinguish two types 
of loops: small ones with winding number w < \/~N, and long ones with w > \/~N. 
The choice w — \/N as boundary is motivated by our observation that the aver- 
age winding number squared (w 2 ), which gives the compressibility of the system, 
makes a jump of the order of N at the condensation temperature. Figure || shows the 
temperature-dependence of the average number (N 00 ) of particles contained in long 
loops, 

JV 

(N°°) = J2 ( N »)> w 
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Figure 3: Fraction of particles contained in long loops. 



as fraction of the total number of particles present. Because of Eq. (j3~l|), revealing that 
the particle number in the string picture is nothing but the total loop length measured 
in units of a = hf3, Fig. |3] can be equivalently interpreted as showing the temperature- 
dependence of the combined length of long loops as fraction of the total loop length. 
We see that long loops are exponentially suppressed in the normal phase, implying that 
here only small loops containing just a few particles are present. The larger the system 
is, the more suppressed the long loops are. Below the condensation temperature, long 
loops appear in the system. We see their number increasing with decreasing tempera- 
ture. At the absolute zero of temperature, virtually all particles are contained in long 
loops. With increasing particle number N, the curves converge to the solid line given 

by 



N 



N- 



(T/T c 



\3/2 



for T <T C 
forT > T c . 



(36) 



implying that also this loop length, like the string tension, exhibits critical behavior 
at T c . The right-hand side of Eq. ( j36| ) is the well-know expression for the number of 
condensed particles. From this we conclude that the particles contained in long loops 
are condensed in the zero-momentum ground state, so that BEC indeed corresponds in 
the string picture to the appearance of long loops, i.e., to the proliferation of strings. 

An interesting question is how many particles are needed for BEC. We investigate 
this point numerically by calculating the heat capacity for various values for the particle 
number N (see Fig. We observe that already for a particle number of the order 10, 
the behavior of the heat capacity at low temperatures closely resembles that of a system 
with many particles. From this we conclude that about 10 particles suffice for BEC to 
arise. 

In the final part of this Letter, we wish to establish contact with Feynman's path- 
integral representation of the partition function of a nonrelativistic system and its inter- 
pretation in terms of rings of cyclicly permuted bosons [§]]. To this end, we recast Eq. 
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Figure 4: Heat capacity for different numbers of particles. 
@ in the form: 

Z = Y,~n E [v Wl+1 G Wl (0)/w 1 ].--[rj w '+ 1 G Wl (0)/w l ], (37) 

i = l Wl ,...,Wi = l 

where the summation over I arises after expanding the exponential function in a Taylor 
series; I physically denotes the number of loops. We omitted the constant term corre- 
sponding to I = 0. It is helpful to introduce the constraint £,- =1 w j = N together with 
an extra summation over TV, so that 

00 OC OC 

E -E E (^) 

Wl ,..., Wl = l N=l wu _. Wl = 1 

in Eq. ([37|). We next group together factors with the same winding number. For given 
I, the right-hand side of Eq. ( |3~7j ) has Z factors, which can be permuted in l\ possible 
ways. If a permutation contains rj identical factors with winding number j, we have 
to divide by rj! because they are indistinguishable. Consequently, 

00 00 oc , . 

z = E 7^ E E TT7T [y^Gm/T • • • [v^a^oy^" 

1=1 L - N=l n U J ] r 

n,...,rjv=0 

£,«■*=' 

= E^i E TT ^[V +1 G 1 (0r-..[^+ 1 G J ,(0)] r *, (39) 
jv=i ' „ 11 J J r 

rx,...,rjsr=0 



N 



E 3 ^ = 

where in the last equation we carried out the sum over I thereby losing the constraint 
Ej rj = I, and inserted unity 1 = N\/N\. The connection with Feynman's path- 
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Figure 5: Interpretation of a loop with winding number w = 3 as the worldlines of 
three cyclicly permuted particles. 



integral representation of Z is established by noting that a propagator with winding 
number j > 1 can be represented as a product of j propagators each of unit winding 
number: 

VGj(0) = J d 3 Xl Jd 3 x 2 --- Jd 3 Xj G 1 [(0,xi) 1 (ft/3 > x a )] (40) 

xGiP,x 2 ), (2/i/3,x 3 )] • ■■G 1 [((j - l)n/3,x,), W^i)], 

where we inserted j ' — 1 intermediate points xa • • ■ Xj, and integrated over them so 
as to allow for all possible values of the intermediate positions. The integral J d 3 xi 
yields the volume factor at the left-hand side of Eq. (|40|). Owing to the periodicity in 
Euclidean time r with period Tif3, Eq. (^5|) may be equivalently written as: 

VG 3 {{))= Jd 3 Xl fd 3 x 2 --- fd 3 Xj Gi[(0,xi),(&Ax 2 )] (41) 
xG x [(0,x 2 ), (?j/3,x 3 )] • • ■Gi[(0,x i ), 



Physically, as was first observed by Montroll and Ward [|13[|, a single loop with winding 
number j is reinterpreted as describing j particles which after a time h(3 are being 
cyclicly permuted: 1 — > 2, 2 — > 3, • • • , j — > 1 (see Fig. ||). 

Amusingly, this picture is reminiscent of Wheeler's attempt to explain why all elec- 
trons have the same mass and charge by supposing that their worldlines form a single 
huge knot [p0|. When the knot is cut by a hyperplane of fixed time, many worldlines 
corresponding to many electrons would appear. The flaw in Wheeler's argument is — as 
was pointed out to him by Feynman — that to each worldline going forward in time there 
should be one going backward in time, corresponding to a positron. This would imply 
an equal number of electrons and positrons, which is not realized in nature. Our picture 
is not nullified by this critique because, owing to the periodicity in r, the worldlines 
discussed here always move forward in time. 

We recognize in Eq. (|9|) apart from the identity permutation G\ (0), an interchange 
of two particles G2(0), a three-particle cyclic permutation Gs(0), and so on. Since 
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Figure 6: Connection between the string (left panel) and cycle representation (right 
panel) of the partition function. 



3 r j = N, the integer N, which was introduced in Eq. (f38|), denotes the total num- 
ber of particles present in the cycles. Figure || gives a two-dimensional illustration of 
how the loop and cycle interpretation are connected. In the left panel, traces of eight 
loops in the i-y-plane are depicted, each of which is parameterized by the Euclidean 
time r. A closer inspection reveals that, for example, the longest loop has winding 
number w = 7. In the right panel, a filled circle is drawn whenever a time %j3 elapsed. 
They mark the starting point of a random walk executed by one particle during the time 
< t < fi/3 and the end point of one executed by another particle. In the longest loop, 
7 particles are cyclicly permuted in this way. 

Now, any permutation of N elements can be factored into cycles of length rjjj = 
1, 2, • • • , N). The number M(r\ , r% , ■ ■ ■ , rjy) of such permutations is given by [|p 

AT! 

M( ri ,r 2r --,r N ) = (42) 

Hi J 3T 3- 



which is precisely the combination appearing in (|39|). In this way, we arrive at the well- 
known path-integral representation of the partition function due to Feynman [|J |J: 

Z= E^i E ^"V) [d 3 xi---d 3 x N G 1 [(0,x a(1) ),(nf3 ) x ai2) )] (43) 

N=l ' aeS(N) 

xGi[(0,x ct(2) ), (7ij9, x <7(3))] • ■ ■G 1 [(0,x a{N) ), (h 13,^(1))], 

where the sum a is over all permutations S(N) of N particles, and 

, 1 .„ , f for bosons 
" =2 (1 ~ ??) = i 1 forfermions. (44) 

(A similar representation of Z, but this time not involving path integrals, was obtained 
earlier by Matsubara pl|].) We thus have established the relation between the string 
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picture and Feynman's, where BEC is connected with the appearance of large rings 
involving the cyclic permutation of many particles. 

In conclusion, we have described a nonrelativistic Bose gas as a path integral over 
spacetime loops. The loops are characterized by their string tension and the number 
of times they wind around the imaginary time axis. The winding number of a given 
loop corresponds to the number of cyclicly permuted particles in a ring in Feynman's 
theory. The string tension, which is determined by the chemical potential, exhibits crit- 
ical behavior at T c . Owing to the vanishing of the string tension, the loops proliferate, 
meaning that for an infinite system, loops with arbitrary large winding numbers appear 
[Bp. Particles contained in long loops were shown to be condensed in the ground state. 
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